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Abstract 

Starting from an array of four-leg Hubbard ladders weakly doped away from half-filling and 
weakly coupled by inter-ladder tunneling, we derive an effective low energy model which contains a 
partially truncated Fermi surface and a well defined Cooperon excitation formed by a bound pair 
of holes. An attractive interaction in the Cooper channel is generated on the Fermi surface through 
virtual scattering into the Cooperon state. Although the model is derived in the weak coupling 
limit of a four-leg ladder array, an examination of exact results on finite clusters for the strong 
coupling t-J model suggests the essential features are also present for a strong coupling Hubbard 
model on a square lattice near half-filling. 

PACS numbers: PACS numbers: 71.10.Pm, 72.80. Sk 
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I. INTRODUCTION 



The microscopic mechanism that generates high temperature superconductivity in the 
cuprates continues to be controversial. One set of proposals is based on the analogy with 
heavy fermion metals where a superconducting dome is observed surrounding the quantum 
critical point (QCP) that arises as antiferromagnetism is suppressed by an external parame- 
ter such as pressure.- In this case the pairing glue arises from the exchange of the soft longitu- 
dinal antiferromagnetic fluctuations in the vicinity of the QCP. In the cuprates doping plays 
the role of the external parameter and there are several proposals for the nature of the QCP 
that appears near optimal doping involving fluctuations in various order parameters e.g. 
nematic,- d-density wave^ and orbital currents^ in addition to antiferromagnetism.- A sec- 
ond set goes back to Anderson's very early proposal that the strong singlet nearest neighbor 
correlations in the 2-dimensional Heisenberg antiferromagnet generates pairing when doped 
holes are introduced. The advocates of this resonant valence bond (RVB) mechanism point 
to the strong asymmetry in the cuprate phase diagram between the physical behavior on the 
under- and overdoped sides of optimal doping and the QCP. This contrasts strongly with the 
symmetric dome observed in heavy fermions. Further the highly anomalous physical proper- 
ties that characterize the pseudogap phase at underdoping are associated with a short range 
spin liquid in the cleanest cuprate materials, e.g. YBa2Cu408 and HgBa2Cu04+2;. Nonethe- 
less strong correlations and the absence of a broken translational symmetry in the pseudogap 
phase have proved to be formidable obstacles to constructing a comprehensive microscopic 
RVB theory for underdoped cuprates. For more details see several recent reviews.-"- 

Several years ago we proposed a 2-dimensional array of weak coupled 2-leg Hubbard 
ladders as an example of a model where occurs a truncation of the full Fermi surface to 
pockets associated with hole or electron doping in a system without broken symmetry.— 
Subsequently this model led to a phenomenological ansatz for the propagator in underdoped 
cuprates starting from a renormalized mean field description of an undoped RVB spin liquid 
insulator.— This phenomenological propagator has been recently used successfully to fit a 
range of experiments covering many anomalous properties of the pseudogap phase.—"— In 
this paper we extend our earlier analysis to the case of an array of lightly doped 4-leg 
Hubbard ladders with an onsite weak interaction. Our goal is to construct a tractable 2- 
dimensional model with a partially truncated Fermi surface in which d-wave pairing arises 
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on the residual Fermi surface through scattering in the Cooper channel. 

Earlier numerical renormalization group studies on the 2-dimensional Hubbard model 
were interpreted as pointing towards a similar pairing mechanism.— A key feature of the 
present model is the presence of a finite energy Cooperon resonance in the pseudogap which 
is generated in association with the partial truncation of the Fermi surface. D-wave pairing 
follows on the remnant Fermi surface through the coupling to the Cooperon. 



II. FOUR LEG HUBBARD LADDERS 

The properties of a single 4-leg Hubbard ladder with open boundaries have been studied 
extensively in both the weak and strong coupling limits. We consider here the former with 
equal nearest neighbor hopping to, along the legs and rungs. In this case the 4 bands split 
into two band pairs. The inner pair, 2, are standing waves on the rungs with wavevectors 
(27r/5, 37r/5). At half filling the corresponding Fermi wavevectors are kp^-^ = ±37r/5 and 
kFA2 = ±27r/5 leading to a common Fermi velocity, vpA = 2to sin(27r/5). The outer band 
pair, -01,2; have Fermi wavevectors Kfb2 = ±7r/5 and KpBi = ±47r/5 and a smaller Fermi 
velocity, vfb = 2to sin(7r/5). 

We obtain a band structure of four bands with energies 

EA.Ak) = e\\{k)T2tocos{2n/5), 

Eb,,M = e||(fc)T2tocos(7r/5), (2.1) 

where e\\{k) represents the dispersion along the ladder. The annihilation (creation) operators 
of electrons of the outer and inner bands, denoted as Bi2,Bl2 and /li2,Al2 respectively, 
are 

4 

Bi = sin(7m/5)c,i; 

71=1 

4 

B2 = sin(47rn/5)cn; 

n=l 
4 

Ai = sin(27rn/5)c„; 

n=l 
4 

A2 = ^sin(37rn/5)cn, (2.2) 

n=l 
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FIG. 1: On the l.h.s. of the figure is picture a four leg ladder with equal hopping along and between 
the legs of the ladder. On the r.h.s are pictured the corresponding four bands, Ai^2, Bi^2 of such a 
ladder. 

where c„ is the corresponding annihilation operator of an electron on the n-th leg of the 
ladder. 

Close to half filling the Fermi velocities of the outer band pair labeled by B are smaller 
than those of the inner bands labeled by A, so that in the presence of interactions the 
effective dimensionless coupling constants for electrons in the inner bands are smaller than 
those for the outer bands. In the weak coupling limit, i.e. an onsite interaction characterized 
by t/ ^ t, this Fermi velocity difference leads to a large difference in the characteristic 
energy scales and to a decoupling of the RG flows of the two band pairs. The outer band 
pair has the larger critical energy scale and flows to strong coupling first as the energy 
scale is lowered.—"— The inner band pair has a lower critical scale. Therefore in the first 
approximation one can treat inner and outer bands of individual 4-leg ladders as decoupled 
from each other. Then each band pair will effectively constitute a two-leg ladder. It is well 
known that two-leg ladders acquire spectral gaps for quite general interaction patterns. For 
the inner bands the smaller dimensionless couplings lead to smaller spectral gaps. At half 
filling each band pair is exactly half filled and behaves as a half filled 2-leg Hubbard ladder. 
The difference in the energy scales leads to a finite doping range x < Xc where all the doped 
holes enter the inner band pair and the outer band pair remains exactly half filled. We note 
in passing that similar behavior is found also in the strong coupling limit, U ^ t.— 

Given that a 4-leg ladder can be reduced to two 2-leg ladders, we will now recall some 
basic facts about 2-leg ladders. For general interactions they become either Luttinger liquids 
or dynamically generate spectral gaps. In the latter case an increased symmetry appears 
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at small energies where a half filled 2- leg ladder can be well described by the 0(8) Gross- 
Neveu model.— The Gross-Neveu model is exactly solvable for all semi-simple symmetry 
groups and a great deal is known about its thermodynamics and correlation functions. In 



the SO (8) case the correlation functions were studied in Refs. (|26 



271 ). Since the model itself 



has Lorentz symmetry, all excitation branches have relativistic dispersion laws: 



E{p) = v^M^ + M^. (2.3) 

The spectrum consists of three octets of particles of mass A and a multiplet of 29 excitons 
with mass y/SA. Two octets consist of quasi-particles of different chirality transforming 
according to the two irreducible spinor representations of SO (8), while the third octet con- 
sists of vector particles. The latter include magnetic excitations as well as the Cooperon 
(a particle with charge ±2e). The 16 kink fields, carrying charge, spin, orbit, and parity 
indices, are direct descendants of the original electron lattice operators on the ladders. 

The SO (8) GN model describes several different phases related to one another by particle- 
hole transformations. Which phase is realized depends on the bare interaction. In this paper 
we assume that it is in the so-called D-Mott phase (in the terminology of Ref. |25| ) . On the 
two two-leg ladders {A and B), the superconducting (SC) order parameters are given by 

Ab = Bi^-^Bi^i — i?2,t-B2,i- (2.4) 

The distinct feature of the half filled ladder is that this order parameter is purely real and has 
a Z2 symmetry. However, the symmetry is restored to U(l) and the phase stiffness becomes 
non-zero as soon as doping is introduced. It is an interesting feature of the SO (8) GN 
model that the only mode which becomes gapless at finite doping is the Cooperon. Neither 
magnetic excitations, nor quasi-particles become gapless.— When the doping increases the 
S0(8) GN model gradually crosses over to the S0(6) GN one plus the U(l) Gaussian model. 
The latter model describes the fiuctuations of the superconducting phase. The effective low 
energy bosonized Lagrangian density for the Cooperon field, $, is 

C=^[v-\djy + v{dJ)% $ = Aoe^t (2.5) 

OTT 



where (j) is the field dual to 6. (Here - according to Ref. (l28l ) - the Luttinger parameter K 
depends weakly on doping and is always in the range 1 > K > 0.9. On the other hand, the 
phase velocity is strongly doping dependent.) 



For values of doping close to the Cooperon band edge — A/2| ^ A) spectral curvature 
is important and the action given in Eqn. 12.51 is inadequate. A better description of the 
Cooperon dynamics is given by the sine-Gordon model 

^ = ^ [vp^dref + VFid^e - 4/x)2] - ^ cos(^). (2.6) 

where = — 4yU^. The mass term here can be thought to arise as follows in a mean 
field way from the SO (8) Gross-Neveu model. The SO (8) Gross-Neveu model can be written 
in terms of fundamental fermions (which are non-local with respect to the original fermions 
in the problem) with an interaction term of the form 

HLT' = Uj:^iryi.^? (2.7) 

a 

Here a = 1, 4 and ipa = ii^ay V'a ) '^^ is a Pauli matrix acting in i? — L space. The four 
fundamental fermions correspond to the different degrees of freedom in SO (8): charge, spin, 
orbital, parity. The Cooperon (charge) we take to be given by tpi. With a finite chemical 
potential lowering the Cooperon gap, the fluctuations of the Cooperon will be strongest. 
Invoking mean field theory, we thus replace i/jlT^iIJa for a = 2,3,4 by its expectation value. 
The resulting bosonization of the remaining degree of freedom ipa results in the sine-Gordon 
model. 



III. SUPERCONDUCTIVITY OF ARRAYS OF FOUR-LEG LADDERS: TWO 
SCENARIOS 

Having elucidated the properties of individual 4-leg ladders, we now consider an array 
of such ladders. We assume initially that the electron-electron interaction acts only inside 
individual ladders and is much smaller than the bandwidth W ~ 2tQ. It is also assumed 
that W ^ tj_ (the inter-ladder tunneling). We imagine two scenarios. In the first we assume 
t± is on the same order as A a, the gap on the inner bands of the four leg ladder, but much 
smaller than A^, the gap on the outer bands. In this case coupling the ladders together 
lead to small Fermi pockets, very much like in Ref.llO). However in this case the pockets 
are found near ±7r/2,±7r/2. The residual coupling between these Fermi pockets and the 
A-cooperons then leads to superconductivity in the A-bands. And because of a proximity 
effect, the superconductivity of the A-bands induces superconductivity in the B-bands. 
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In the second scenario, we assume A a ^ tj. ^ A^. In this case t± wipes out the 
effects of interactions on the A-bands. Coupling them together then gives us an anisotropic 
two dimensional Fermi liquid. But as is much smaller than A^, the Cooperons on the 
outer bands at zeroth order remain unperturbed. The coupling then between the anisotropic 
Fermi liquid and the B-Cooperons induces superconductivity in the system as a whole. This 
superconductivity is d-wave in nature. 

We now elaborate on these two scenarios. 



A. Scenario I 



We treat the interladder hopping through a random phase approximation (RPA) analysis 
of the interladder hopping. The form of the hopping is taken to be long range 

-f^interladder = ~ ^a]b ^n,a^rn,b, (3-1) 

nytm,a,b 

where a,b = 1, ...,4 run over the legs of an individual ladder and n and m mark the n'th 
and m'th ladders. By particle-hole symmetry the hopping is assumed to have peaks both 
near k± = and k = G/2 where G = (0,7r/2) is the inverse lattice vector perpendicular to 
the ladders. In particular the hopping takes the form 

t:T = i^-i-^r-nfab{m-n) (3.2) 

where fabijn — n) = fba{n — m) and fabiS^) = (i.e. no (additional) hopping within a ladder). 

By treating -ffmteriadder in an RPA approach, we find that the single particle Green's 
function takes the form 



RPA(a;, k,, ky) = V^{ky)Gf^{u, k,, ky) + V^{ky)Gf^{U, k,, k 



Gf{u,k,,ky) = GaX^^^^J (3.3) 

' l + GAX^,k^W^{kyy 

where 

t\ff{ky) = 2j2^0Si4ky) 

n>0 

2(Si + S2)tij + (2SiS2 - Si)(ti;2 + t2,l) - 2SiS2(t3;i + ^1,3) - 52(^1,4 + kl ) 1 
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n-1 n n+1 

123412 3 412 3 4 



n+1, n'**:^- ' ' . ' ' 

t4j=i,...,4 

FIG. 2: An array of four leg ladders. As an example of the hopping assumed in the RPA analysis 
(hopping between every second leg), we show how electrons can hop between the fourth chain of 
the n-th 4-leg ladder and the chains on the n + 1-th ladders. 



€i^{k,) = 2 J]cos(4A:,) 

n>0 




Ti{ky) = 2{sl + si) + 2{2siS2 — si) cos{ky) — 4:SiS2Cos{2ky) — 2slcos{3ky); 

r2iky) = 2{sl + si) - 2(2siS2 - si) cos{ky) - 4siS2 cos{2ky) + 2sl cos{3ky), (3.4) 

and Si = sm(7r/5) and S2 = sin(27r/5). We have assumed the hopping is real and that 
the low energy contribution to G^^ comes from the A-bands as A a ^ A^. Thus 
G Ai{oo , kx) / G A2{^ , kx) are the Green's functions of the A-band electrons on a given 4- 
leg ladder. As we have discussed in the previous section Gai/Ga2 are no more than the 
bonding/anti-bonding electron Greens functions for a 2-leg ladder. The RPA does not mix 
Gai and Ga2 as the weights of the two are found near differing Fermi wavevectors (i.e. we 
can take GAi{k)GA2{k) ~ safely for all k). The presence of ri{ky) and r2{ky) act as 
structure factors which cause the quasiparticle weight at various ky to be negligible. While 
the denominator of G^^ has the periodicity of the reduced Brillouin zone i.e. ky and 
ky + 7r/2 are identified) these structure functions merely have the periodicity of the original 
zone i.e. ky and ky + 27t are identified). 
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FIG. 3: The electron and hole pockets of an array of weakly coupled four leg ladders shown in a 
periodic zone scheme. On the I.h.s. of the figure are pockets at zero chemical potential. On the 
r.h.s. of the figure are pictured the pockets for finite chemical potential such that the interladder 
hopping satisfies 2Aa + 2fi > \to\ > 2Aa — 2fj.. 

The Green's functions for A1/A2 at zero chemical potential are given by 

where the E^^ are defined in Eqn. 12.11 At a chemical potential, /x, that does not exceed the 
gap, Ga, is given by GaX^^ /^^ k) = Ga,{uj - /i, 0, k) 

The excitations are then given by the locations of the poles in G^^ ^P-^^ These poles then 
imply that the excitations have the dispersion relation 



Ei{k^,ky) = /i - 



± J {EAAK)-tf^ iky) /2y + A\ 



(3.6) 



For sufficiently large t^^-^ a Fermi surface forms (found by solving Ei = 0) consisting of 
electron and hole pockets. The type of pocket is determined by the sign of the effective 
hopping 

tl^^ {ky) > 2Aa + 2/i — )■ electron pocket; 
tf^iK) < -2Aa + 2/i -> hole pocket. (3.7) 

In our conventions a positive chemical potential favors hole pocket formation while disfa- 
voring electron pockets. As t1^^ [ky) grows beyond this minimal value, the pockets grow in 
size. We take the hopping such that 



t/\ky Ky) 



-to(l - {ky - Kyf/nl + ...), Ky ~ 0, ±7r/2 
to(l - {ky - Kyf/nl + ...), Ky^ ±7r/4, ±37r/4 



(3.8) 



where kq is the small parameter guaranteeing that the RPA is a good approximation. 
The dispersion relations of the quasi-particles near the hole pockets are 

E,{k,, ky) = ^^l^Zlil. + - ^y^' _ e^. (3.9) 

where p. = ±K,^^ T e,. = 7. = (1 - ^(A^ " /^^ + K"mY^ 

fT^±i = eff, s , and m\\i = ' ' -2 — -. 

In Figure 3 are plotted the expected Fermi pockets. On the l.h.s. of Figure 3 are plotted 
the pockets found at zero chemical potential while on the r.h.s. are plotted the pockets for 
a chemical potential such that 2Aa + 2fi > to > 2A^ — 2yU. For such a condition one obtains 
only hole pockets. We see that hole pockets occur in the vicinity of (±7r/2, ±7r/2). 



1. Luttinger sum rule 

The Luttinger sum rule (LSR) for the single particle Green's functions at the particle 
hole symmetric point takes the form 

2 



n 



i2n) 



d'^k, (3.10) 

Giuj=0,k)>0 



where n is the electron density. The corresponding Luttinger surface of G{uj, k) is defined as 
the loci of points in /c-space where G{u) = 0, k) changes sign. These sign changes occur both 
at the poles and the zeros of G. In order to apply the Luttinger sum rule, we must take 
G{uj, k) to be one of G"^^^, i.e. we must apply the LSR to each band separately (see Eqn. 
(13. 3p for the definition of G'jf ). (We only apply the LSR to the electrons in the A-bands 

1/2 

- the LSR also holds separately for electrons in the B-bands.) 

At the particle-hole symmetric point, zeros are present in G^^^_^(0,k) along the lines 
ky = ±Kf^^. In the absence of pockets the LSR is satisfied G{u = 0,kx = KF^^^,ky) 
because of these zeros. And when t*^-^^ becomes strong enough so that pockets form, the 
appearance of equally size electron and hole pockets on either side of ^Kp^- ensure that 
the Luttinger sum rule continues to hold. 

Introducing a finite chemical potential (with fi < A a/ 2) leaves the LSR violated as 
expressed in Eqn. fl3.10p . However it continues to hold in a modified form. Because in a finite 
chemical potential, the ladder Greens functions are given by GAi {u}, /i, k) = GAi — /i, 0, k), 
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the LSR holds if we consider the sign changes the Green's function undergoes not at a; = 
but at a; = /i. 



2. Superconducting Instability 

The residual interactions between the Fermi pockets and the Cooperons will lead to 
instabilities in the RPA solutions as temperature goes to zero. Provided a finite chemical 
potential is present the leading instability will be to a superconducting state. While gapless 
quasi-particles only exist in the A-bands, both A and B bands will go superconducting 
simultaneously. The general form of the Cooperon-quasiparticle interaction is 



1 ^(q,k,k') ^ 

9 NTn ^QPAi^^^)^QPA{^.<^) 



2 ^ NLa 

q,k,k' 

A^p^(k,q) = e..,[AL(k + q)<,(-k)-AL(k + q)AL,(-k)]. (3.11) 

Here L is the length of the ladders, a is the interladder spacing, and is the number of 
ladders in the array. ^a,b are the Cooperon fields whose bare propagators are defined as 

DK..,k) = (rt,(k..„)tt(k,.„)>„ = (3.12) 

We see that g has the dimensionality of energy x length^ and Fj has the dimensionality of 
energy x length^/^ . 

The different terms in Eqn. (13.111) have different origins. The strongest interactions are 
presumably Fj as this term already exists for uncoupled ladders. Inter-ladder interactions, 
such as interladder Coulomb repulsion, also contribute to Fj. However interladder hoping 
does not - this contribution is suppressed due to a mismatch between the Fermi momenta 
of the Ai and Bi bands. The coupling g is smaller than Fj: it arises only in second or- 
der perturbation theory from intraladder interactions and from presumed weak inter-ladder 
Coulomb interactions. 

The pair susceptibility for the quasiparticles Aqpa in an RPA approximation is given by 
Xgpi(^n,q) = jj^Y. rrfre--(TAQp^(ki,q,r)A^p^(k2,q,0)) 



ki,k2 
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2C{uJn, q) 



(3.13) 



1 + g{q)C{u^, q) - 2 r2(q)C(a;„, , q) q) ' 

We have assumed that the couphngs 5f(q, k, k') and rj(k, q) are such that we can ignore 



their dependence on k and k'. Here C(a;„, q) is the Cooper bubble: 



C{uJn, q) = 2 



dJvx dky 
47r2 



/(e^i(k + q))-/(-e^i(-k)) 



(3.14) 



eAi(k + q) -eAi(-k) 

Here eA^/al^) the bare dispersions of the A1/2 quasi-particles. As T — 0, C{un,q = 0) 
develops a logarithmic divergence: C{uJn,q = 0) ~ y^mjjrnX log (^^^). 
The pair susceptibility for the Cooperons fields has a similar RPA form: 



Xr^(a;„,q) = (T0,(g, r)0t(g, Q)) = A°K,^) + (A°K, fc))2r2(g)xg^^(a;„, q) 

DO(a;„, k) + g{q)C{u, q) - 2D?(a;„, q)C(a;„, g)r?(g) 



(3.15) 



1 + (?(q)C(a;„, q) - 2 r2(q)C(a;„, , q) q) 

where A = 5, i? = A. The superconducting instability occurs when the denominator in 
Eqns. (13.31 and 13.51) vanishes at cj„ = 0, g = 0, that is 



c(o,o)(^7(o)-2^r,^(o) 



0. 



(3.16) 



We note that this vanishing occurs simultaneously in all channels, li g > (though inter- 
ladder Coulomb repulsion is repulsive, the interactions between quasi-particles on a given 
ladder is attractive leaving the sign of g indeterminate) the instability occurs only when the 
chemical potential approaches sufficiently close to A^/2 so that the resulting effective inter- 
action becomes attractive. This chemical potential corresponds to minimal doping at which 
the superconductivity appears. Taking ^ A^, the corresponding transition temperature 
takes the form 



Tc = max{Tci, Tc2}; ^ epi exp 



1 



2ri(o)^_ 

(Ai-4/.2) 



(3.17) 



where epi 



4 ^P(0)-2m 

pockets, the density of dopants is equal to 



If we suppose that fj, and t*^-^^ are such that we only have hole 



x(/i) 



(2A, + |tf^(0)| + 2/i)(-2A, + |t,^ 



eff 



2/i) 



27/2^^2 1 1^// (0)1 1/2 



tf^(0)|+2/x 



(3.18) 
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FIG. 4: The Fermi surface of the A-bands in a periodic zone scheme. 

If we denote the critical doping as xdfi = where the A-Cooperon becomes soft, we 

see that the transition temperature behaves as ~ exp(— a(xc — x)) as x approaches Xc, 
that is to say, the transition temperature has a strong dependence on doping. It should 
be emphasized that this critical doping Xc as defined above does not coincide with the 
optimal doping as typically understood. Optimal doping can be thought of as the doping 
level associated with a change in the Fermi surface topology. However in this understanding 
our model always remains in the underdoped regime since the quasiparticle Fermi surfaces 
remain small as far as the interladder tunneling remains much smaller than the gap of the 
outer (B) band pair. 

B. Scenario 2 

We now consider the second scenario where As ^ t± ^ A a- Because t± is much larger 
than A A but smaller than A^, the effects of the interactions are wiped out in the A-bands 
while preserved in the B-bands. In particular a gapful Cooperon still exists on the B-bands 
while the coupled A-bands appear as an anisotropic two dimensional Fermi liquid. 

We can distinguish two parameter ranges in this scenario. At small dopings fi < Ab/2, 
the B-Cooperons remain gapped. The effective Hamiltonian for the two dimensional Fermi 
liquid in the A-bands and the Cooperons in the B-bands appears as 

if^^ = 5^ei(k)4(k)A(k) + 62(k)4(k)A2(k) + 5^i?^^(A;.)4(k)$^(k); 

k k 

6,(k) = EAXh)+tf\ky) 
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EbM = ^Jkl + Al-2fi, (3.19) 

where EaX^x) is given in Eqn. (12. ip and fj^{ky) in Eqn. (13. 4p . We illustrate the two 
dimensional Fermi surface of the A-bands in Figure 4. 

The form of the quasi-particle-Cooperon interaction is that of Eqn. (I3.13P (though of 
course, now we have no A-Cooperon and so this coupling is absent). This system, like in 
Scenario 1, has a pairing instability to superconductivity. The pairing susceptibilities in an 
RPA approximation take a similar form as for Scenario 1: 

X ^ 2C(a;„, q) 

XqpaK^u, qj ^ ^ ^(q)C(a;„, q) - 2r|(q)C(a;„, , q)D° (a;„, q) ' 

^RPA(,^ „N _ D'l{uJn.k)+g{q)C{uj,q) 

^ "'^^ l + ^(q)C(u;„,q)-2r|(q)C(a;„,,q)Z}U^n,q)- ^ ^ 

where C(a;„,g) is defined as in Eqn. (I3.15p . 

As we no longer have pockets as in Scenario 1, but instead have an anisotropic 2D Fermi 
liquid whose Fermi surface consists of slightly deformed lines (see Figure 4), the divergent 
with temperature behaviour of C(0, 0) now takes the form 

^^(O,0) = ^log(^^^^). (3.21) 

Because the A-quasi-particles are already gapless, a finite /i dopes the A-bands with doping 
x'^(/i). Thus EfaXi^) = EfaXi^ = 0) ~ fi^vpA- If we denote the critical doping, /ic, as 
the doping when the B-Cooperon becomes soft (i.e. fic = Ab/2) and x-^ = x^{fic) the 
corresponding doping of the A-bands, we can rewrite the form of the B-Cooperon propagator, 
DO (0,0), as 

D^iun = 0, g = 0) ~ , ^. (3.22) 

Again we emphasize that the critical doping x^ as defined above does not coincide with 
optimal doping - in this model we are always in the underdoped regime. For this range of 
doping we obtain a transition temperature of the form 



rr I Nl/2 / 2r|(0)fi;'B fr.^ 

Tc = {eFAieFA2) ' exp - -navFAX „.2 „iii^.a\i ~ 9W 



(3.23) 



Vpj^a'^{{x^y — x^ 

and we see that the critical temperature grows extremely fast with doping, similar to the 
transition temperature determined in Scenario I. 
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The second region occurs at x > Xc, when the holes penetrate into the outer B-bands. 
Here the 0(8) Gross- Neveu model governing the B-bands undergoes a crossover into a 0(6) x 
f/(l) Gross- Neveu model. The B-Cooperon propagator aX u,k = becomes more singular. 
At the same time the velocity of the phase fluctuations becomes small and these fluctuations 
can be treated as slow modes. Integrating over the nodal fermions one obtains the effective 
Lagrangian for the phase fluctuations: 



£ = ^ [- Jccos Q(0„(a;) - 0„+i(x)) j 



+ ^ (vF{f^)idJn - 4/i)^ + VF{fl){dren)'^ - y COs(^), (3.24) 

where n is a sum over ladders. As we have already noted the parameter K is renormalized 
by the Coulomb interaction to be slightly less than 1. vplfi) is more dramatically affected, 
taking the form vf{^) ~ ^'Fsi-^ — 1)^^^ so that it vanishes at x = Xc (or equivalently 
fi = Ab/2). As a side remark we note that there is an alternative way of presenting the 
effective Hamiltonian. The above Lagrangian (Eqn. I3.24p is the continuum limit of the 
following model: 

H — I ^ •J{'''n,m+l'^n,m + h.C.) — Jc{T^+l^rn'^n,m + h.C.) + 

n,m 

[(-l)"M-2/i]<„}, (3.25) 

where are Pauli matrix operators. In the continuum limit r~ becomes the order parameter 
field e*2. Here J ~ M. The model presented above is a model of anisotropic spin-1/2 
magnet on a 2D lattice with a staggered (M) and uniform magnetic fields (2/i). This form 
of the Hamiltonian has been proven to be very convenient for numerical calculations yielding 
promising results for the transport.— 

We again estimate the transition temperature using an RPA argument. At T = the 
doping of the entire system (both the A and the B bands) is 

a: = /xpA + c^(|^-l)^/^ (3.26) 



where c is a constant and oa = — - — The detailed 

avpATT • — 



a single chain at T=0 can be extracted from Ref 



30 



"orm of the Cooperon propagator for 
However to obtain an estimate for 



Tc, it is enough to use the finite temperature Luttinger liquid expression for the Cooperon 
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propagator: 



1/2K ^ 



DI(t,x) = z( ^— ^ — , , / , , --, (3.27) 

where Z is a numerical constant. Thus 

DlicU = 0, q = 0) ~ _ iy/2^1/2Kj.-2+l/2K^ (g ^g) 

Substituting the latter expression into RPA expressions for the pairing susceptibilities (Eqns. 
I3.20p we obtain an estimate for the critical temperature upon doping as follows: 

A 2-2K 

T.~(x-^)^^-\ (3.29) 

This dependence on the doping is much weaker than (Eqn. I3.27p . It holds in the region 
where phase fluctuations are already strong. 

Thus we have obtained two regimes with different doping dependence of T^. The first 
one is the BCS-like with given by Eqn. fl3.23p . It corresponds to lowest doping levels. 
The other regime, which in our model still describes a situation an anisotropic ID-like Fermi 
surface, is the regime with strong phase fluctuations. The mean field transition temperature 
in this regime is given by Eqn. fl3.29p . A further increase of doping presumably will lead to 
a change in the Fermi surface topology and is not considered in this paper. 

IV. DISCUSSION 

Phenomenological models based on coupled fermions and bosons similar to that derived 



here, have been proposed much earlier Refs.(l32H35l) to describe the high temperature su- 



perconductors. The closest similarity are to the models proposed by Geshkenbein, loffe and 
Larkin^ and by Chubukov and Tsvelik.— Both these phenomenological models examined 
Fermi arcs centered on the nodal directions, coupled in the d-wave channel to Cooperons 



associated with the antinodal regions. The model studied in Ref. (1341 ) had dispersionless 
Cooperons which provided BCS-style coupling for the nodal quasiparticles. The result was a 
superconducting transition with weak fluctuations, similar to our x < Xc case. In the model 



considered in Ref. 



35l the Cooperons possessed a one-dimensional dispersion which resulted 

34 -36 con- 



in strong fluctuations as takes place in our case for x > Xc- The authors of Refs 
sidered the fluctuation regime above T^ when the Cooperon energy is close to the chemical 
potential and drew comparisons to experiments in several underdoped cuprates. The key 
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ingredients controlling superconductivity in the array of 4-leg Hubbard ladders that we have 
considered in this letter, are a small residual Fermi surface (either pockets as in Scenario I 
or arcs as in Scenario II), which is coupled in the d-wave Cooper channel to a finite energy 
Cooperon associated with the pseudogap responsible for the partial truncation of the Fermi 
surface. The properties of a weak coupling 4-leg Hubbard ladder near to half-filling are 
used to obtain these key ingredients. Our goal is to derive a tractable model containing the 
important features that are relevant to high temperature superconductivity in the cuprates. 
In order to assess the relevance of our model to this goal, clearly one must examine whether 
these key ingredients are present in a two dimensional Hubbard model on a square lattice 
near half-filling. 

As we mentioned above, earlier numerical renormalization group studies on the 2- 
dimensional Hubbard model were interpreted as pointing towards a similar pairing mecha- 
nism arising from enhanced pairing correlations present in a condensate that truncates the 
Fermi surface in the antinodal regions. There are of course two reservations in these earlier 
works. Firstly, the one loop approximation in the numerical renormalization group stud- 
ies limits them to at most moderately strong onsite repulsive interactions. Secondly, the 
renormalization group studies per se break down when the scattering vertices flow to strong 
coupling and the nature of the resulting low energy or low temperature effective action is 
a difficult problem which could only be surmised rather than explicitly derived. These two 
weaknesses make it imperative to examine the question whether these key ingredients are 
present also for strong coupling. 

The most reliable strong coupling calculations are exact diagonalization studies of strong 
coupling Hamiltonians. The only limitation is the flnite cluster size which currently is limited 
to small clusters containing up to 32 sites and 1,2 and 4 holes. Leung and his collaborators 
Refs. 



37 



39 



have reported a series of calculations for these clusters using the strong coupling 
t-J model and its extensions to include longer range hopping and interactions. We begin 
with a recap of the main conclusions of these calculations. The allowed set of k-points in a 
32-site cluster with periodic boundary conditions contain both the four nodal (7r/2, 7r/2) and 
two antinodal points (7r,0) and (0,7r). A single hole enters at a nodal point. For two holes 
there are two different states that are possible groundstates depending on the parameter 
values. For the plain t-J model with only nearest neighbour hopping a 2-hole bound pair 
state with d(x^ — y^) symmetry is the groundstate on the 32-site cluster for J/t > 0.28. The 
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binding energy is quite small at J/t = 0.3 but grows with increasing J/t. An extrapolation 
from finite size clusters to the infinite lattice however suggests that the pair state is no 
longer the groundstate at J/t =0.3, but an excited state with an energy of approximately 
0.17t.— The inclusion of longer range interactions and hopping in the t-J model increases 
the energy of the pair state further and confirms the conclusion that for parameter values 
relevant to cuprates the groundstate of the cluster has two unbound holes in the nodal 
states.— Extending the calculations to the 32-site clusters with 4 holes, which corresponds 
to a doping of 1/8, shows all 4 holes entering into nodal states with no signs of pairing 
correlations.— In view of the prominent bound pair excited state for 2 holes, a low energy 
excited state with two of the holes in a bound state may also be expected here. However at 
present there is no information on this question to the best of our knowledge. 

Leung and collaborators^!"^ concluded from these calculations that at low densities holes 
entered the nodal regions, possibly in pockets, and as a result there was no evidence for d- 
wave pairing correlations in the groundstate for realistic values of the parameters in t-J 
models. However the analysis presented here suggests a more optimistic conclusion. First 
we note that the nodal points in the 32-site cluster are very special, because exactly at these 
points the coupling in a Cooper channel to a d-wave Cooperon vanishes by symmetry. Thus if 
we interpret the d-wave pair excited state as evidence for a finite energy Cooperon in the t-J 
model and its extensions, then as the occupied holes at finite doping move out from the exact 
nodal points, a d-wave pairing attraction is generated through the coupling to this Cooperon, 
similar to the scenarios we discussed earlier. Note an earlier study for two holes on smaller 
clusters by Poilblanc and collaborators^ concluded in favor of the interpretation of the 2-hole 
bound state as a quasiparticle with charge 2e and spin 0, which would be an actual carrier 
of charge under an applied electric field. In other words they concluded that a Cooperon is 
present in the strong coupling t-J model at low doping. A more detailed analysis of the origin 
of the pairing in this state was published recently by Maier et al.— Note the hole density in 
the case of 2 holes in a 32-site cluster is very low so that the superconducting order we are 
postulating should coexist with long range antiferromagnetic order. There is considerable 
evidence both numerical, in variational Monte Carlo calculations, and experimental, in favor 
of such coexistence, as discussed in the recent review by Ogata and Fukuyama.- 

We conclude that there is strong evidence that the pairing mechanism in the present 
model is not confined to weak coupling and ladder lattices, but will also operate in the 
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strong coupling t-J model on a square lattice at low doping. 
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